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Exploring degeneracies in modified gravity with weak lensing 
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By considering linear-order departures from general relativity, we compute a novel expression for 
the weak lensing convergence power spectrnm nnder alternative theories of gravity. This comprises 
an integral over a ‘kernel’ of general relativistic quantities multiplied by a theory-dependent ‘source’ 
term. The clear separation between theory-independent and -dependent terms allows for an explicit 
understanding of each physical effect introduced by altering the theory of gravity. We take advantage 
of this to explore the degeneracies between gravitational parameters in weak lensing observations. 


I. INTRODUCTION 

In recent years, weak gravitational lensing has been put 
forth as a promisingmethod of testing gravitation on cos¬ 
mological scales [iHa] , with some exciting first constraints 
having been found already Moreover, advances in 

relevant data analysis (for example, [HI) and the coming 
next generation of lensing-optimised surveys mean that 
we will soon be in a position to take full advantage of the 
potential of weak lensing. 

Stronger constraints on gravity are obtained by com¬ 
bining weak gravitational lensing with other probes. One 
observable which is commonly touted as providing par¬ 
ticularly complementary constraints to weak lensing is 
fasia). Here, /(a) is the linear growth rate of structure, 
defined as: 


This paper is structured as follows: in Section im we 
detail the derivation of the expression for Pk(£). Section 
m discusses how weak lensing degeneracy directions be¬ 
tween gravitational parameters can be understood with 
the help of our expression. Finally, Section HVl provides 
forecast constraints on deviations from GR-I-ACDM from 
future surveys, and interprets these constraints using our 
expression for Pk.{£)- We conclude in Section IVl 

II. CONVERGENCE IN MODIFIED GRAVITY: 
THE LINEAR RESPONSE APPROACH 

In what follows, we use the scalar perturbed Friedmann- 
Robertson-Walker (FRW) metric in the conformal New¬ 
tonian gauge, with the following form: 


/(a) 


din Am (a) 
din a 


( 1 ) 


where Am (a) is the amplitude of the growing mode of 
the matter density perturbation, and erg (a) is the am¬ 
plitude of the matter power spectrum within spheres of 
radius 8 Mpe/h. The combination /erg (a) can be con¬ 
strained through measurements of redshift-space distor¬ 
tions in g alaxy surveys. 

In [Q, an expression for fas{a) was derived in the 
case of linear deviations from the model of general rela¬ 
tivity (GR) with AGDM. Here we build on this work by 
constructing a similar expression for Pk(-^), the angular 
power spectrum of the weak lensing observable conver¬ 
gence (k). The main advantage of our expression is that 
it clearly distinguishes the physical source of all modified 
gravity effects to Pnii), which allows for a more thor¬ 
ough interpretation and understanding of these effects 
than previously possible. While we focus on Pk.{£) in 
this work, recall that the two main weak lensing observ¬ 
ables, convergence and shear, can be trivially intercon- 
verted [l^. Therefore, we treat convergence as a proxy 
for weak lensing more generally, and all expressions which 
we derive could be equivalently and easily formulated in 
terms of shear. 
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ds^ = a{T)'^ [—(1 -I- 2'I')c?t^ + (1 ~ 2^)dxidx'^] . (2) 


Our parameterisation of alternative theories of gravity 
makes use of the quasistatic approximation (see, for ex¬ 
ample, H)- The quasistatic approximation states that 
within the range of scales relevant for current galaxy sur¬ 
veys, the most significant effects of a sizeable class of 
modihed theories can be captured by introducing two 
functions of time and scale into the linearised field equa¬ 
tions of GR. These functions play the role of a modified 
gravitational constant, and a non-unity (late-time) ratio 
of the two scalar gravitational potentials: 


2V^$(a,fc) 
$(a, k) 
\I'(a, k) 


8TTGa^fi{a, k)pMAM{a, k) 
l{a,k). 


(3) 


In GR, both 7(0, k) and /i(a, k) are equal to I. 

Glearly equation [3] can only be an effective description 
of more complicated, exact sets of field equations [H- 
[^ . However, several works have numerically verified the 
validity of the quasistatic approximation in many gravity 
theories (notably those with one new degree of freedom) 
on the distance scales considered here [251 - 1^ . 

We first compute the power spectrum of the conver¬ 
gence in general relativity, and then generalise to al¬ 
ternative theories of gravity. We make the simplify¬ 
ing assumption that radiation can be neglected for all 
redshifts of interest in this paper. That is, we take 

= 1. 
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A. Calculating convergence: general relativity 


B. Calculating convergence: modified gravity 


The convergence, k, describes the magnification of an 
image due to leasing. This effect is captured by the 
geodesic equation for the displacement of a photon trans¬ 
verse to the line of sight. In the cosmological weak leasing 
context of general relativity, this is given by: 

(4) 

where , b indicates a partial derivative with respect to 0^, 
X is the radial comoving distance, and 
is a two-component vector representing on-sky position. 
This equation can be integrated to obtain the ‘true’ on- 
sky position of the light source as a function of the ob¬ 
served on-sky position. The convergence is then given by 
taking the two-dimensional on-sky Laplacian (V^) of this 
expression: 

KGRiO) = x)9{x) (5) 

where g{x) is the leasing kernel: 

g{x) = 2 x 1^°° dx' (l - ^ (X'), ( 6 ) 

W (x) is the normalised redshift distribution of the source 
galaxies, and Xoo is the comoving distance at a —t 0 . 

We compute the power spectrum of the convergence 
following closely the method laid out in . In the small 
angle approximation, it is straightforward to find: 

i ir 

^ / (0 (7) 


where x labels three-dimensional position such that x = 
{xdi,xd2,x) and x' = ( 0 , 0 , xO- * and j label the source 
redshift bins to be considered. 

Performing the integrals over 0i and 02 and then over 
ki and ^ 2 , we have: 

K’^i)=I dx ^ dx' gAx') 

^ JO X Jo 

“ ( 8 ) 

Finally, the Limber approximation [Uls^, valid here on 
Z ^ 10 [I|, is employed, such that A :3 <C and therefore 
k Ki -k. The small angle limit also means that Pl;^{l) — 
where i labels an angular multipole [1^. We find: 





5*(x)Sj(x) □ 



(9) 


We have computed here the power spectrum of the con¬ 
vergence; that of the shear could be straightforwardly 
calculated by replacing equation [S] with the appropriate, 
similar definition. 


As indicated in equation [3l generally in non-GR theo¬ 
ries $ 7 ^ fk. So, in modified gravity equation 0] becomes: 


— (X^') =-(<!>,f,+ «'.&). 

The convergence then becomes: 
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( 10 ) 


Rmg = 


1 


dx V" 


dx 


'H{x) 


$(6»,x)+ ^'(d,x) gix) 


$( 6 »,a:)-h 4 '( 6 »,a:) g{x{x)) 


( 11 ) 


where hereafter we will use x = ln(a) instead of x or 
a, and we have converted the integration measure to x 
using dx = —c/'Hdx, where "H = aH is the conformal 
Hubble factor. Note that x here is distinct from the 
three-dimensional position variable x. 

To calculate the power spectrum of the convergence un¬ 
der modifications to GR, we follow [l^ and perturb our 
field equations about those of the GR-I-AGDM model. 
Our reasoning here is that current observations only per¬ 
mit theories which can match GR-I-ACDM predictions to 
leading order; we are interested in determining next-to- 
leading order corrections that are still permitted. Note 
that we are building a theory of linear perturbations in 
model space, which is distinct from spacetime perturba¬ 
tion theory. We define the perturbations of the qua¬ 
sistatic functions /i and 7 about their GR values using: 

pi{x, k) = \ + 5pL{x, k) 

'){x,k) = \-\-5x(x,k). ( 12 ) 


In addition, we introduce a perturbation about the stan¬ 
dard value of the effective equation of state of the non¬ 
matter sector, w{x): 

w{x) =-I-\-(i{x), (13) 


and we define the useful related quantity: 

u{x) = f fJ{x')dx'. (14) 

Jo 

We now consider how these linear perturbation vari¬ 
ables propagate through to k and hence Firstly, 

from equation [3l we can write: 

<l)(x, k) + '^[x,k) = (1-1- , ^ ) $(a;, k) 

\ l{x,k)J 

Ki {2 — Sx{x,k)) ^{x,k). (15) 

In order to express our results as corrections to 
GR-I-ACDM, we need to relate ^{x,k) to ^GR{x,k). 
There are two effects to be accounted for. Firstly, the 
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relationship between $(a;, k) and matter density pertur¬ 
bations can be altered. Secondly, if the field equations 
are modified, Amixyk) will evolve at a different rate, 
and hence will be displaced from its GR value. To ac¬ 
count for this we introduce the deviation 5^^(x,k) = 
AM{x,k)/A%^{x,k) - 1. In [13 it was shown that 
Sa{x, k) is given by the following integral expression: 


where we have defined Gi{xix)) = ■ 

There are still two non-GR effects to account for, 
both originating from the modified expansion history. If 
/3(x) ^ 0 in equation [CT 'H(x) and x(x) will scale differ¬ 
ently with the time variable x. Using the expression for 
6'H{x) ='H{x) — 'Hgr{x) derived in equation IA4I in the 
Appendix, we find that 


(5a (a:, k) = 


il^{x)I{x, x) dSf{x, k) dx. 


(16) 


The integrand above separates into two parts: 6Sf{x,k), 
which encapsulates all deviations from GR-I-ACDM, and 
il^{x)I{x, x), which is a weighting function containing 
GR-I-ACDM quantities only. It will be useful for us to 
present the explicit form of SSf{x, k) here, derived in [l^ : 

dSf{x, k) = 6n{x, k) — Sj{x, k) 


+ 


(1 - 


3 nff (1 -b fanix)) u{x) + fGR{x)/3{x) 


(17) 

The explicit form of I{x,x) can be found in M- 

With these modifications in hand, the parameterised 
Poisson equation becomes: 

—2fc^$(a;, k) = p%P{x)Am{x, k) (1 -b 5p,[x, k)) 

= 3HGfi(2^)^M”(2;)AM^(a;, k) 

X (1 -b (5A(a:, fc)) (1-b i5/i(a;, fc)) (18) 

where in going from the first to the second line, we have 
used the fact that the combination 'Hg,^(a;)D£/^(x) is un¬ 
changed by our modifications to the background expan¬ 
sion rate, as shown in Appendix Hence, $(x, k) is 
given in terms of ^gr{x, k) by: 

<I>(a:, k) ~ ^gr{x, k) (1 -b k) + Sp,{x, k)). (19) 

Combining equations [15] and HOj we now have an ex¬ 
pression for <I> -b 'k in modified gravity in terms of the GR 
potential plus perturbative correction factors: 


$(0:, k) + 'I'(a::, k) ~ ^gr{x, k) ^2 — 5^{x, k) 
-b 25 a( a;, k) + 26p{x, k)j . 
So, referring to equation ITTl k becomes: 

1 


( 20 ) 


Rmg{(^) = 


4^ 

X (2 




■H{x) 

2Sfi{x, k) — S"f{x, k) 


-2SAix,k)) . ( 21 ) 


At this stage, it becomes more convenient to work di¬ 
rectly with the power spectrum (£). This can be com¬ 
puted to linear order in deviations from GR-bACDM, in 
direct analogy to the method outlined for the GR case in 
Section Hi Al We find: 

^4 rO 


p:’^(£) = 


4 

X (1- 


dx- 


'H^x) 

■ 25p,[x, k) — 5x{x, k) 


Gi(x{.x))Gj{x{x))Pi^ 

2Sa{x, k)) 


1 


nix) 


nGRix) 

1 


1 - 


5n 


nGRix) 


nGRix) 

and hence 


1 - -j^uix)il - n%^ix)) 


(23) 


Xix) 

5xix) 


nGRix') 

fO 


1 - 


snix') 


dx' 


nGRix') 


nGRix') , 

■pruix') (1 - ^^ix')) dx' (24) 


where 5% = x - XGR- 

The deviation of x(a^) from its GR value will also affect 
quantities which depend on xix)^ such as Gixix)) and 
P§^i^/xix)) [13■ We allow for this by expanding these 
in a Taylor series around XGR, to first order: 


P 


GR 


Xmg 


P 


GR 


XGR 


91nP<i. 


ainx 


XGR 


<5x 

XGR 


G^ix)GAx) 


GiiXGR)GjixGR) 

( ainG*(x) 


1 -b 


- 


91nx 


d\^G,ix) 

91nx 


(25) 


<5x 

^ XGR 
' (26)“ 


where 5x is given by equation above. We have now 
accounted for all modified gravity effects, and these are 
summarised in Table |T| 

Finally, it will be more convenient for us to work in 
terms of Ps, the matter power spectrum, instead of P^. 
We do so via the following expression, where for clarity 
we temporarily omit the label ‘GR’ on all quantities.: 




9.\[ix)Dix)'^Ps ix = 0,k). 


(27) 


Here Pix) is the usual growth factor of matter pertur¬ 
bations. Inserting equation 1271 into equation 1251 we find: 


P 


GR 


xix 

( 22 ) 


,xix) 


Pi 


1 - 


GR 


XGR 

d ln(fc“^ P^^ix = 0,k)) 


dink 


^=i/XGR 


5x 

XGR 


( 28 ) 








































4 


Correction 


Description Equation 


^{x,k) + ^{x,k) ~ (1 + ::^^)${x,k) 

^{x, k) ~ ^gr{x, fc)(l + 5A(a:, k) + 5/i(a;, k)) 

x(x) ~ xgr(x) + lf rj^^^ulx) (1 - n%^{x)) dx 


Gi{x)Gj{x) 


Gi{XGR)Gj{XGR) 


( 81nGi(x) I 81nGj(x) \ 
ainx ainx ) 


Non-unity ratio of scalar potentials [15] 

Altered Poisson equation 1191 

Altered 'H{x) [23] 

Altered x [H 

Altered G(x) [26] 




L 

XGR 



IdGR ( £ 

'] pGR ( e '] y; 


Sx 

Altered 

M 

{x 

\xgr)^ 

dink 

XGR 

k=^/XGR 


TABLE I: Here we summarise the various corrections to the GR expression for including a brief description and the 

number of the equation in which they are introduced. 


Drawing together, then, equations [22l [2311261 and [28l and convergence power spectrum under modifications to gen¬ 
using equation 1271 we obtain our final expression for the eral relativity: 

I 




,2 


X 


+ 


fdluGYx) 

, dluGYx) 

d\n{PY^{x = 0, fc)/fc'^)\ 

Sx{x) 

V i^lnx 

dlnx 

91nfc ) 

Xgr{x) 

XGR 2 


1 -h -u{x) (1 - Vt^{x)) 

(29) 


The major advantage of equation [S^ is that it neatly sep¬ 
arates the convergence power spectrum into the familiar 
GR expression (the non-bracketed quantity) and a cor¬ 
rection factor (the bracketed terms). It is then easy to 
pick out contributions from: 

• the modified clustering properties (described by d/i 
and dy). 


Here we have defined the ‘kernel’ term: 


]C(x t) = ^ 9i{XGR{x))gj{XGR{x) 


16 


^ D'gr{x) 


Xgr(.xY 
^Git(^) O' 


AGfl(a:) 


GR( 
M \ 


(31) 


• the modified expansion history (described by /3, u and the ‘source’ term: 
and (5x), and 


• the modified growth rate of matter density pertur¬ 
bations (encapsulated in (5 a, see equation [T6l) . 

It will be useful for us to write equation [29] in a form 
which explicitly highlights the GR expression and the 
correction factor: 

P^'J(£)= / dx IC{x,£){l + SSwLix,£)). (30) 

J —OO 


5Swl{x,£) = ^u{x) (1 - ^M^{x)) 

+ 2Sfi(x, k) — Sj(x, k) + 2 Sa{x, k) -b 


( dlnGiix) 
y ainx 


ainG^ 

91nx 


dln[PY^{x = 0,fc)/fc^] 

dink 


XGR 


Sx{x) 

Xgr{x)' 

(32) 
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III. UNDERSTANDING DEGENERACIES straightforward algebra allows us to find an expression of 

WITH THE LINEAR RESPONSE APPROACH the form 


We have at hand an expression (equation [29]) for 
under modifications to general relativity. Let us now in¬ 
vestigate what this can teach us about the degeneracies 
between gravitational parameters in weak leasing obser¬ 
vations. Note that we restrict ourselves to discussing de¬ 
generacies between parameters describing modifications 
to gravity. We do not examine degeneracies between 
gravitational and cosmological parameters, nor do we in¬ 
vestigate degeneracies with the galaxy bias. We leave 
these questions for future work. 

In this section, we consider the case in which and 
S'y are independent of scale, due to the fact that the 
scale-dependence of these functions is expected to be 
sub-dominant to their time-dependence [Tj, [l^, ■ We 

will briefly investigate scale-dependence later, in Section 
II V Cl Additionally, as we are working in the quasistatic 
approximation, our analysis is restricted to the regime of 
validity of linear cosmological perturbation theory. Var¬ 
ious values of imax which ensure this to be true are sug¬ 
gested in the literature (see for example fli). Adopting 
a conservative approach, we select imax = 100 here and 
for the remainder of this work. 

We first remind the reader of how degeneracy direc¬ 
tions may be calculated. Then, using equations [T7| and 
1291 we explore how the degeneracy directions of weak 
lensing and redshift-space distortions in the space of the 
parameters of 5^{x) and <57(0;) are affected by the cho¬ 
sen ansatzes for these functions. Note that here and for 
the remainder of this work, we compute the GR-I-ACDM 
matter power spectrum using the publicly available code 
GAME [3^ and using the best-fit ACDM parameters of 
the 2013 Planck release (including Planck lensing data) 

M- 


A. Calculating degeneracy directions 


Degeneracies exist when an observation can probe only 
some combination of the parameters we wish to con¬ 
strain. The degeneracy direction is the relationship 
between parameters in the fiducial scenario (here, in 
GR-l-ACDM). For example, if this relationship is a = b, 
then the relevant observation can probe only a — b, not 
a or 6 individually. 

In the case of weak lensing, degeneracy directions can 
be understood in the following schematic way. First, de¬ 
fine the fractional difference between (£) in an alter¬ 
native gravity theory and in GR-I-ACDM: 




pyr.M) 


(33) 


To hnd the degeneracy direction, we find the relationship 
which exists between parameters when SP^’^{i) = 0. If 
we consider a two parameter case (call them a and 6), 


a = V{l)b (34) 

where 'D(P) may be a complicated expression, but de¬ 
pends only on GR-I-ACDM quantities. The degeneracy 
direction, we see, depends on i in the weak lensing case. 

In order to calculate 6Pl’^{i), we need to specify Wi, 
the normalised redshift distribution of the source galaxies 
in the redshift bin i. We select a source number density 
with the following form: 

n(z) oc z“e (^0) , (35) 

and we select a = 2, P = 1.5, and zq = 2;m/1.412 where 
Zm = 0.9 is the median redshift of the survey, mimicking 
the number density of a Dark Energy Task Force 4 type 
survey Si- In this section, we will simply consider all 
galaxies between 2; = 0.5 and z = 2.0 to be in a single 
redshift bin, with W (y) given by normalising equation 

ESI 

To break parameter degeneracy in a two parameter 
case, a second observable with a different (ideally orthog¬ 
onal) degeneracy direction is introduced. Here, we choose 
this second observable to be redshift-space distortions, as 
it is known to provide nearly orthogonal constraints to 
weak lensing. We will therefore often employ results from 
[1^ . Particularly, we reproduce here their equation for 
the deviation of fcrs{x) from its GR value, analogous to 
our equation ESI 

6fas[x) = - -=/ Gf{x,x)6Sf[x)dx 

J ^8 \^/ J — 00 

(36) 

where 6Sf{x) is given as in our equation]!?] and Gf(x, x) 
is a general relativistic kernel given in equation 34 of [l2j| . 
Note that fa^{x) above is independent of fc, because we 
are considering a case where /r and 7 are functions of 
time only. 

The degeneracy direction of a measurement of fa^{x) 
can then be computed in a directly analogous way to that 
described above for weak lensing. The sole difference is 
that instead of depending on multipole £, the degeneracy 
direction is dependent on the time of observation, x. 

With this information in hand, we now explore de¬ 
generacy directions of weak lensing and redshift-space 
distortions in the space of the parameters of 6fi{x) and 
Sj{x). 


B. Degeneracy directions in the po — So plane 

As mentioned above, redshift-space distortions are the 
preferred choice of an additional observation to break 
weak lensing degeneracy in this scenario. Upon closer ex¬ 
amination, this statement hinges upon the chosen time- 
dependent ansatz for the functions which parameterise 
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deviations from GR. As there is no clear front-runner 
amongst alternative theories of gravity, typically a phe¬ 
nomenological ansatz is chosen, in which deviations from 
GR become manifest at late times in order to mimic ac¬ 
celerated expansion. It is for this type of phenomenologi¬ 
cal ansatz that redshift-space distortion and weak leasing 
observations are known to provide complimentary con¬ 
straints 0. 

However, it may also be desirable to constrain the pa¬ 
rameters of a specific theory of gravity. The functions 
which parameterise the deviation of an alternative grav¬ 
ity theory from GR can, in principle, take on a wide range 
of time-dependencies. Is the combination of weak leas¬ 
ing and redshift-space distortions still an effective way to 
break degeneracies and constrain the parameters of the 
theory we consider? A priori, this is unknown. 

To explore this issue, we consider now the degeneracy 
directions of weak leasing and redshift-space distortions 
under two different ansatzes for the functions which pa¬ 
rameterise deviations from GR. For this section only, we 
make the simplifying assumption that /3(a:) = 0 (i.e. the 
expansion history is ACDM-like). We expect that the 
effect of this assumption on our qualitative findings will 
be small. 

First, we perform a simple operation on 6^{x) and 
< 57 ( 0 ;) to obtain a more observationally-motivated set of 
functions. Let us call these Ji{x) and S(x), in keeping 
insofar as possible with the notation used in Q. The 
choice of this set of functions allows nearly orthogonal 
constraints in the /2o — Sq plane for the phenomenologi¬ 
cal choice of time-dependence. The mapping between the 
two sets of functions, as shown in Appendix |B1 is given 
by: 

i;(a:) = 5^i{x) - i(57(x) 

fjL{x) = Sfj,{x) — S"f{x). (37) 

We can rewrite the linear response ‘source’ terms for both 
weak leasing (equation [32]) and redshift-space distortions 
(equation |3n|) in terms of fl{x) and S(a;) (in the I3{x) = 0 
case): 

/ X 

Qfp'{x)I{x, x) Jl{x)dx 

-00 

5Sf{x) = fi{x). (38) 

We see that SSf{x) depends solely on fl{x). 

The expression for 5Swl{x) requires slightly more 
pause. It depends on E(a;), but it also contains another 
term, which comprises an integral over p.{x) and some 
general relativistic quantities. By comparing with equa- 
tion[T6l we can easily recognise this term as 26 This 
term quantifies a correction to the degeneracy direction 
of weak lensing away from Eq = 0. It is clearly depen¬ 
dent upon the ansatz of time-dependence chosen for p,{x). 
Particularly, we note that due to the integral nature of 
the correction term, choices of p.{x) which persist signif¬ 
icantly over longer times will result in greater deviations 
to the degeneracy direction. 


We now consider two ansatzes for fi{x) and E(x). First, 
consider a phenomenological ansatz, for which we know 
weak lensing and redshift-space distortions to be an effec¬ 
tive combination in constraining gravity theories. This 
choice is a specific case of the form proposed in and 
has been used in, for example, It is given by: 


fi{x) = flo 


E(a:) = Eq 


nf^{x = 0 ) 
= 0 ) 


(39) 


where O^^(x) is the time-dependent energy density of 
dark energy in the fiducial ACDM cosmology. 

We insert 6Swl{x) (equation|35|) into equations I551and 
|36|with our chosen fl{x) and E(a;). We then follow the 
procedure sketched in Section llll Al to find the degeneracy 
directions of weak lensing and redshift-space distortion in 
the fiQ — Eq plane. In this particular case the degeneracy 
direction of redshift-space distortions does not depend 
on time. This is because SSf{x) is dependent on only 
one parameter, p,o, and therefore the only degeneracy 
direction is /to = 0. 

The degeneracy directions for this ansatz can be seen 
in Figure ID (left). In the case of weak lensing, we have 
plotted the degeneracy direction for £ = 50; directions for 
other multipoles £ = 10 — 100 differ only within 5%. We 
see that, indeed, the degeneracy directions are nearly or¬ 
thogonal, with only a slight correction of the weak lensing 
degeneracy direction away from Eq = 0. 

Now, consider selecting an ansatz with a very different 
time-dependence. To guide our selection, recall that we 
expect choices of p,{x) which persist over longer times to 
result in a greater value of the integral term in equation 
1551 and hence a greater deviation of the weak lensing 
degeneracy direction from Eq = 0. Therefore, with no 
attempt to correspond to any particular gravity theory, 
we select the simplest possible choice which persists over 
long times: constant /I(x) and E(a:). 


E(a:) = Eo 

p.{x) = p. 0 . (40) 


In reality, we use step functions beginning at 2 ; = 15 
rather than true constants to allow for the numerical 
computation of the degeneracy directions. The degener¬ 
acy directions are calculated as before, and are plotted in 
Figure |T](right). Clearly, they are less orthogonal than in 
the previous case, as expected from the comments above. 

What does this example tell us about the effective¬ 
ness of combining weak lensing and redshift-space dis¬ 
tortions? The ansatz for p,(x) and S(x) given by equa¬ 
tion sni deviates from GR-I-ACDM at all times after 
2 ; = 15. As mentioned above, most cosmologically- 
motivated alternative theories of gravity present devi¬ 
ations from GR-I-ACDM at late times only, mimicking 
accelerated expansion. Therefore, we treat the case of 
equation 1401 as a heuristic ‘upper bound’ on the cumula¬ 
tive effect produced by the integral term of equation 1381 
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The effect of this term can be quantified by considering 
the angle of the weak leasing degeneracy direction with 
respect to the vertical. We find that for the range of ^ 
which we consider and for the ansatz given by equation 
|40l the maximum possible value of this angle \s 9 ~ 50°. 
Although the degeneracy directions in this case are cer¬ 
tainly no longer orthogonal {9 = 0°), they are sufficiently 
distinct that we expect the resulting constraints to be 
reasonable (if not ideal). We have therefore shown that 
the effectiveness of combining weak leasing with redshift- 
space distortions in the /3(x) = 0 case is relatively robust 
to the chosen form of /2(x) and 


IV. FORECAST CONSTRAINTS FROM 
FUTURE SURVEYS 


In addition to providing an understanding of degenera¬ 
cies, our expression for enables the forecasting of 

constraints. The straightforward form of eauation l30l ren- 
ders the calculation of Fisher matrices very simple, and 
clarifies the interpretation of the resulting forecasts. We 
take advantage of these features to forecast constraints on 
gravitational parameters for a Dark Energy Task Force 
4 (DETF4) type survey, as defined in the classification 
of [15 . We focus on combined constraints from weak 
lensing and redshift-space distortions, with some consid¬ 
eration given as well to baryon acoustic oscillations. 

As mentioned above, the forecasts presented here em¬ 
ploy the technique of Fisher forecasting (see, for example, 
[39|)- The key quantity of this method is the Fisher in¬ 
formation matrix: 


Pab = — 


lnC \ 

dpadpb / 


(41) 


where pi are the relevant parameters, and C is the like¬ 
lihood. For redshift-space distortions, we straightfor¬ 
wardly build on the results of [ll] to construct the appro¬ 
priate Fisher matrix. However, for weak lensing we re¬ 
quire a slightly different expression. Although our equa¬ 
tion allows for the cross-correlation of source galaxy 
redshift bins, we have until now considered only a sin¬ 
gle wide redshift bin. In practice, weak lensing data are 
normally considered in a number of tomographic redshift 
bins. In (dOj, the Fisher matrix for such a situation is 
shown to be given by: 


Pab — ^ 




/sfeyTr [Cqj^C^, 


PgrCA 


(42) 


where ,a is a derivative with respect to Pa, fsky is re¬ 
lated to the fraction of the sky observed {fsky = 0.375 
for a DETF4-type survey), and C is an Nf, x Nb matrix 
where Nt is the number of tomographic redshift bins. C 
represents the observed power spectrum of the conver¬ 
gence, and is given by the following expression [^ : 

{£) = (£) + (43) 

Ui 


where fii is the number density of galaxies per steradian 
in bin i and is the rms intrinsic shear, equal to 0.22 

for a DETF4-type survey. 

Computing the appropriate value of fii requires the 
selection of source redshift bins. In practice, once in pos¬ 
session of data, the selected bins are those which are max¬ 
imal in number while maintaining shot noise sufficiently 
below the signal. For our forecasting purposes, we in¬ 
stead follow, for example, Q and Q- We select redshift 
bins by subdividing n{z) of equation [35] into 5 sectors, 
such that the number of galaxies in each bin is equal. 
The value of fi for the total redshift range for a DETF4- 
type survey is given by h = 3.55 x 10®, so the value in 
each tomographic bin is simply fii = h/5. 

In the following subsections, we use the Fisher formal¬ 
ism to compute forecast constraints in a number of sce¬ 
narios. We first consider constraints on the parameters 
of 6fi{x) and 6j(x) in the case where we fix the expansion 
history to mimic ACDM. We then incorporate expected 
measurements of wq and Wa from Baryon Acoustic Oscil¬ 
lations to forecast constraints on the parameters of Sp{x) 
and d"f{x) in the case where we marginalise over the pa¬ 
rameters of /3(x). We finish by discussing the directions 
of best constraint in the parameter space of the scale de¬ 
pendent ansatz for p,{x, k) and ^{x, k) put forth in Cl. 


A. ACDM-like expansion history: /3(a;) = 0 


We first consider constraints on the parameters of dfi{x) 
and Sj(x) in the case where the expansion history is fixed 
to be ACDM-like. As in equation|37l we transform 6pL{x) 
and Sj(x) to p,{x) and S(x), and we choose the time- 
dependence given by equation 1331 

To compute these constraints, we calculate the 2x2 
Fisher matrix for lensing, for redshift-space distortions, 
and for both observations combined. For this, we require 
expressions for the derivatives 


d 

dpo 

d 


p::^{£) + 


fii ) 


Pl;^{£) + 


fii J 


dpo 

5So 


(44) 

(45) 


These are found in a straightforward manner from equa¬ 
tions |29| and [361 we present them in Appendix jC] 

The resulting forecast constraints are illustrated in Fig¬ 
ure |2l As discussed in Section Hill the degeneracy direc¬ 
tions of the two observables are nearly orthogonal in this 
case. Combining them results in promising forecast con¬ 
straints on po and Eq. We see that we can expect a 
DETF4-type survey to provide constraints at a level of 
approximately 4% in this plane, in the case where f3{x) 
is assumed to be fixed at 0. 
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FIG. 1: Degeneracy directions of weak lensing {£ = 50, dashed red) and redshift-space distortions (solid green) in the no — So 
plane, where /r(a;) and E(a:) scale as Da^(®) (left)i s^nd as constants (right). 



FIG. 2: Forecast constraints for weak lensing (orange), 
redshift-space distortions (green) and both observables com¬ 
bined (blue) for a DETF4-type survey, in the po — So plane 
with I3{x) fixed to 0. Gontours represeirt the 68.3% and 95.4% 
confidence regions. 


B. The effect of marginalising over {uio,uia} 


In reality, /3(x) is not fixed to zero, but rather the as¬ 
sociated parameters will also be constrained with some 
non-zero error. While weak lensing and redshift-space 
distortions will provide some constraints on these, it is 
baryon acoustic oscillation (BAO) measurements which 


are expected to provide the best constraints on the ex¬ 
pansion history of the universe. 

In this section we use a CPL-type ansatz for /3{x) as 
proposed in [U [4^: f]{x) = wq + 1 + Wa{l — e^). We 
incorporate forecast BAO constraints on wq and Wa and 
use these to obtain expected constraints in the fig — So 
plane. We hrst marginalise over only wg, while holding 
Wa to its fiducial value of 0; then we examine the effect 
of allowing Wa to vary as well. 


1. Marginalising over wo; Wa = 0 

We first demonstrate how constraints in the fto — ^o plane 
are affected by marginalising over wg when Wa is held 
fixed to its hducial value of 0. 

Because we are now incorporating information about 
three parameters (/tq, Sq and wg), our Fisher matrices 
are 3 x 3 in dimension. In order to compute these, we 
now require additional derivatives of and fa%{x) 

with respect to wg\ all are listed in Appendix [Cl Because 
transverse measurements of BAO are independent of non¬ 
background gravitational effects [i^, the Fisher matrix 
of BAO is non-zero only in the (wg^ wg) component. The 
value of this matrix component is equal to ^ ^ -, where 

'^■wq.BAO 

Wo,BAO is the 1-cr error on wg from BAO measurements. 

To explore the effect of marginalising over wg, we con¬ 
sider three levels of constraint from BAO; 

1. For comparison: the case where wg is fixed to its 
fiducial value. This is identical to the case consid¬ 
ered in Section HV Al 

2. The case where (Jwo,bao = l9o- This scenario mim¬ 
ics best-case constraints from a DETF4-type sur¬ 
vey. 
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3. The case where Cwo.bao = 5%. This lies between 
current best constraints and scenario 2 above. 

The resulting constraints from the combination of weak 
leasing, BAO, and redshift-space distortions are shown 
in Figure|3]and FigurelH Figure [3] shows in the left-hand 
panel the forecast constraints on /Iq for cases 1 — 3 above 
when marginalising over wo and Sq ; the right-hand panel 
displays the same for Sq when marginalising over wq and 
fto- Figure ID shows the 68.3% forecast joint constraints 
on flo — So in scenarios 1 — 3 while marginalising over wo 
only. 

We note from Figure 0] that the degeneracy direction 
of the combined constraint in the p,o — Eg plane changes 
considerably between the three scenarios. and Eg are 
mildly negatively correlated in scenario 1, whereas in sce¬ 
nario 2 they are positively correlated, and in scenario 3 
even more so. This can be understood by considering 
the joint forecast constraints in the /fg — wq and Eg — wq 
planes, marginalised in each case over the other non-Wa 
parameter. These are displayed at a 68.3% level in Fig¬ 
ure [S] for scenario 3. Both /2g and Eg are shown therein 
to exhibit a positive correlation with wq . This implies 
that p,Q and Eg are also positively correlated with each 
other, except in the case where wq is fixed or constrained 
so tightly that this effect is negated. As the constraint on 
Wq is loosened, moving from scenario 1 through scenario 
2 to scenario 3, this positive correlation becomes more 
pronounced. 

We notice also from Figure|3]that the constraint on Eg 
is relatively insensitive to the level of BAO constraint on 
wg, whereas the constraint on /ig changes considerably 
between scenarios 1 — 3. This is consistent with Figure 
01 in which we see that the degeneracy direction in the 
/2g — Wo plane has a far greater positive slope than that 
in the Eg — wq plane. These degeneracy directions, and 
hence the relative sensitivity of p,o and Eg constraints 
to Wo constraints, can be understood by considering the 
expressions for (equation 011) and 8fa^{x) (equa¬ 
tions [T71 and 011) . Both and Sfas{x) are given by 

integrals in time over a kernel and a source term. In the 
case of Sf(7g,{x), the general relativistic kernel Gf{x,x) 
is significant back to ^ ~ 15, whereas in the weak leasing 
case, the kernel is non-zero only as far back in redshift as 
the furthest source galaxies (2 = 2 in this case). In the 
current model of /3(x), deviations from a ACDM expan¬ 
sion history are more significant at early times, whereas 
p-{x) and E(a;) are both chosen to be significant only at 
late times (below z ~ 5). Therefore, the Sfas{x) integra¬ 
tion from z ~ 15 favours sensitivity to the background 
expansion variable wo over flo , whereas the weak leasing 
integral, significant only from z ~ 2, results in relatively 
greater sensitivity to Eg. This results in the relative sen¬ 
sitivity of the flo constraint to the wo constraint level, 
as seen in Figure 0] Note that we have not accounted 
here for any uncertainty in galaxy bias models at high 
redshifts, which may have significant effects on the sen¬ 
sitivity of fcTsix) to the background expansion at early 
times. 


Finally, we note that there is clearly a direction in the 
/Ig —Eg plane which is entirely insensitive to the change in 
Wo- This is in fact expected due to the nature of the con¬ 
tours displayed. Given the hypothetical 3D confidence 
region in the space of /Ig, Eg and wg, the marginalised 
constraint of scenario 3 is equivalent to projecting this el¬ 
lipsoid into the flo — Eg plane. When we reduce the error 
in only the wo direction as in scenario 2 - that is, re¬ 
ducing the error in the direction orthogonal to the plane 
of projection - the resulting projection will, by simple 
geometrical considerations, coincide with the first pro¬ 
jection in two locations. The same argument can then 
be extended to the case of fixed wo, which involves sim¬ 
ply taking a slice of the 3D ellipsoid at the location of 
the flo — Eg plane. 

2. Marginalising over {wo,Wa} 

We now consider the case where we do not fix Wa to 
zero. In this scenario, there is information present about 
4 parameters (^g. Eg, wg, Wq), so all Fisher matrices are 
4x4. In addition to the previous derivative expressions, 
we now need derivatives with respect to Wa of Pf.'^ (£) and 
fas{x). Once again, these are computed from equations 
051 and 051 and listed in Appendix ICl In this scenario, the 
BAO Fisher matrix is slightly more complicated, as the 
entire 2x2 block related to wo and Wa is non-zero. 

In analogy to the above, we consider three scenarios: 

1. The scenario where wq and Wa are fixed to their 
fiducial values. Again, this for comparison, and is 
identical to the case considered in Section IIV A1 

2. The scenario where the BAO Fisher matrix rep¬ 
resents the best-case expected constraints from a 
DETF4-type survey. In this scenario, the compo¬ 
nents of the BAO-only covariance matrix (the in¬ 
verse of the Fisher matrix) are given by: Cwq^wO = 
0.0010, C^a,wo = -0.0038, and C^a,wa = 0.016 
0 . 

3. The scenario where the BAO-only covariance ma¬ 
trix is obtained by multiplying the covariance ma¬ 
trix listed above in scenario 2 by an overall factor of 
(8.2)^. This corresponds to the case where the pro¬ 
jected 68.3% error on wq from BAO is 5% and all 
other elements of the covariance matrix are scaled 
up accordingly. 

The left-hand panel of Figure 0] presents the combined 
weak leasing, redshift-space distortion and BAO forecast 
constraints on /2g while marginalising over wq, Wa and 
Eg; the right-hand panel does the same for constraints 
on Eg while marginalising over wo, Wa and flo- Figure [71 
meanwhile, presents the 68.3% confidence regions in the 
flo — Eg plane while marginalising over wq and Wa- 

We see that the forecast constraint on Eg is now 
slightly more sensitive to the level of BAO constraint 
on Wo and Wa than in the above case where Wa is fixed. 
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FIG. 3: Forecast constraints from weak leasing, redshift-space distortions, and BAO in the case where wo has been marginalised 
over and Wa has been fixed to 0. The left-hand panel shows the confidence region for flo when Eq is marginalised over, while 
the right-hand panel shows the confidence region for Eq with /to marginalised over. Black, solid: wo fixed; red, dashed: BAO 
error on wo = 1% (DETF4); green, dotted: BAO error on wo = 5%. 



So 


FIG. 4: Forecast 68.3% confidence regions in the fio — Bq 
plane, marginalising over wo, for the case where Wa = 0. 
Black, solid: wo fixed; red, dashed: BAO error on wo = 1% 
(DETF4); green, dotted: BAO error on wo = 5%. 


This is particularly noticeable in scenario 3, in which the 
expansion history is the least well-constrained. Turning 
to /2o, we see from Figure IH that the forecast constraint 
remains sensitive to our knowledge of the expansion his¬ 
tory in much the same way as in the Wa fixed case. That 
is, the constraint in scenario 3 is broadened consider¬ 
ably relative to that in scenario 2, and both are slightly 
broader than in the above case where Wa fixed. Finally, 
examining the combined plot in Figure [71 we see that the 
confidence regions therein are slightly larger than those 
in the corresponding Figured where Wa is fixed (other 


than for scenario 1, which is the same in both figures by 
design). 

We surmise that allowing for a time-dependence in 
the equation of state of the effective dark energy com¬ 
ponent (via /3(x)) loosens the expected constraints on /2o 
and So, but not catastrophically so. In fact, the level 
of constraint provided by BAO measurements on the ex¬ 
pansion history of the universe appears to have a greater 
effect on forecast constraints in the fig — Sq plane than 
does our assumption regarding the time-dependence of 
that expansion history. 


C. Scale-dependent and 'y{x,k) 

Until this point, we have neglected any scale-dependence 
of /i(a:, k) and 7(1, fc), focusing only on time-dependence. 
We now consider a scale-dependent ansatz. 

It has been shown that in the quasistatic regime and 
for local theories of gravity, k) and 7(0;, k) can be 
expressed as a ratio of polynomials in k with a specific 
form [3: 


l{x,k) 

ti{x,k) 


pi{x) +p 2 ix)k'^ 

1 +P3{x)k‘^ 

1 +P3ix)k‘^ 

P4{x) +P5{x)k'^ ' 


(46) 


This form has recently been considered in [3|j in which 
a Principle Component Analysis was undertaken for a 
combined future data set including weak leasing and 
galaxy count measurements from the Large Synoptic Sur¬ 
vey Telescope (LSST), as well as Planck measurements 
and upcoming supernova data. Therein, the primary goal 
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FIG. 5: Forecast constraints at a 68.3% level from weak leas¬ 
ing, redshift-space distortions and BAO measurements, in the 
case where o-wo,bao = 5%. The top panel displays these 
joint constraints in the flo — wo plane, marginalising over Eo, 
while the bottom panel does the same in the Eo — wq plane, 
marginalising over /tq. In both cases, Wa is fixed to 0. 


was to provide insight into future constraints on /r(x, k) 
and ^{x, k) as well as on the overall scale-dependence. 
Here, we use the Fisher matrix formalism to find the level 
of forecast constraint on various parameter space direc¬ 
tions from a DETF4-type survey. Our expression for 
will then allow us to understand the related degeneracy 
structure of the parameter space. 

Because we are interested only in broadly understand¬ 
ing the best-constrained directions of the parameter 
space, we fix /3(x) = 0 for simplicity. In GR+ACDM, 
pi{x) = pa{x) = I and P 2 {x) = p^{x) = Pb{x) = 0. 
Therefore, in keeping with our linear response approach, 
we define: 

pi{x) = 1-1- 5pi{x) 

P2{x) = 5p2{x) 

Psix) = Spsix) 

Piix) = 1 + Spiix) 

Psix) = Sp5(x). (47) 

Then, by dropping higher-order terms in dpiix), we find 


expressions for pix, k) and 7(0;, k): 

p{x, fc) ~ 1 + 5psix)k‘^ — Spiix) — 5psix)k‘^ 

7(x, fc) ~ 1 + Spiix) + 6p2ix)k'^ - 5p3ix)k'^. (48) 

Our previous choices for Spix, k) and S^ix, k) have been 
consistent with eouation ldSl we have simply assumed that 
Sp 2 ix), Spsix)^ and Spsix) have been sub-dominant, and 
let —Sp4ix) = Spix) and Spiix) = Sjix). 

It will be useful to have the equivalent expressions for 
pix, k) and S(a:, k). We find them using equation 1571 

pix, k) ~ —5piix) — dpiix) 

+ i-Sp2ix) + 2Sp3ix) - Spsix)) 

T,ix,k) ~ --6piix) -Sp4ix) 

+ [^\5p2ix) + ‘^P3ix) - Spsix)^ . (49) 


As usual, in order to forecast constraints using Fisher 
matrices, we must select an ansatz for the functions 
dpiix). From equation 1461 we see that in fact, while 
6 piix) and 6 p 4 ix) are dimensionless, Sp 2 ix), 6 p 3 ix) and 
dpsix) must have dimensions of length squared. We 
choose a form similar to that of equation [39] for both sets 
of functions, introducing a mass scale where necessary to 
account for their different dimensionalities: 


^P(iA)ix) =Po 


(1.4) 


^’K^ix) 


Q^^ix = 0) 


SP{2,3,5)ix) = 




n^^ix) 


(20iTo)2 Hf«(x = 0) ■ 


(50) 


Thus, the parameter space is five-dimensional, with pa¬ 
rameters {pI,Pq,Pq,Pq,Pq}. The scale c/20Ho has been 
chosen as it provides sensible numerical eigenvalues. We 
will see below that the numerical value of this scale is 
irrelevant to our results. 

To compute the Fisher matrices we require derivatives 
of and fas with respect to the parameters Pq. These 
are computed in a straightforward manner from equa¬ 
tions |25| and |3ni and are presented in Appendix |C] As 
can be seen there, or as is obvious from equation 1461 
some of these derivative expressions are dependent upon 
k. In the weak leasing case, this is trivially dealt with, as 
the Limber approximation allows us to set k = How¬ 
ever, the derivatives of fas with respect to Pq, Pq and Pq 
are truly /c-dependent. 

In order to treat this case, we abandon the usual (gen¬ 
eral relativistic) notion that all scale-dependence of /erg 
is factored out, and instead divide our forecast observa¬ 
tions into five bins in fc. As in [H , we let these bins 
have edges k = [0.005,0.02,0.05,0.08,0.12,0.15], stop¬ 
ping short of entering the regime where nonlinearities 
dominate. We select the error on the measurement in 
each A:-bin based on the assumption that a DETF4-type 
survey, which covers a large fraction of the sky, will re¬ 
sult in tighter measurements of large-scale modes than 
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FIG. 6: Forecast constraints from weak lensing, redshift-space distortions, and BAO in the case where wo and Wa have been 
marginalised over. The left-hand panel shows the conhdence region for fto when Eq is marginalised over, while the right-hand 
panel shows the confidence region for Eq with p.o marginalised over. Black, solid: wo and Wa fixed; red, dashed: scenario 2 
described in Section TlV B 21 green, dotted: scenario 3 described in Section TlV B 21 



So 


FIG. 7: Forecast 68.3% confidence regions in the fio — Eo 
plane, marginalising over wo and Wa- Black, solid: wo and 
Wa fixed; red, dashed: scenario 2 described in Section TlV B 21 
green, dotted: scenario 3 described in Section TlV B 21 


of small-scale modes. To model this, we divide the total 
error budget of each redshift bin, given by [1^, into the 
hve /c-bins listed above. The first redshift bin (between 
k = 0.005 and k = 0.02) receives 4% of the error budget, 
the next two bins {k = 0.02 — 0.05 and k = 0.05 — 0.08) 
receive 12% of the error budget each, and the final two 
bins (fc = 0.08 — 0.12 and k = 0.12 — 0.15) receive 36% 
of the error budget each. 

We can now construct the Fisher matrices for weak 
lensing and redshift-space distortions, recalling in the 
redshift-space distortions case to sum over bins of scale as 


well as of redshift. We invert the total combined Fisher 
matrix to obtain the covariance matrix, and diagonalise 
to obtain five eigenvectors with corresponding eigenval¬ 
ues. 

Before examining and interpreting these eigenvalues 
and eigenvectors, we pause to consider more fully the 
implications of the difference in dimensionality of 6 pi (x) 
and 6 p 4 {x) vs Sp 2 ix), Sp 3 {x) and 5p^{x). In equation 1501 
P2, Ps and pg have essentially been scaled by 
This scaling factor is arbitrary, and using a different 
scaling factor would alter the the numerical value of 
the eigenvalues of the covariance matrix. Therefore, al¬ 
though all eigenvalues of the covariance matrix are di¬ 
mensionless, any forecast constraint on P2i Pa or P%^ 
or on any combination thereof, must be multiplied by 
(20go)^ appropriate scaling factor). This provides 

a physically meaningful value with dimensions length 
squared, which represents a constraint on the combi¬ 
nation P(2 3 5) ^ ( 20 HoY • Without introducing a prior, 
it is impossible to meaningfully compare numerical con¬ 
straints on parameter combinations within two different 
parameter sub-spaces: that of p\ and p® and that of p®. 
Pa: and pg. 

First, consider those eigenvectors in the subspace of p® 
and P4. We list them here, along with the square root of 
the associated eigenvalue This acts as a measure 

of how well the parameter space is constrained in the 
direction of the eigenvector. We have: 

cfi = 0.8pi — 0.004p2 + O.OOlpa -I- 0.6p4 -I- 0.002p5 

~ 0.8pi + 0.6p4 (51) 



















13 


with ~ 0.01, and 

di = —0.6pi + 0.0003p2 + 0.0002p3 + 0.8p4 — O.OOOGps 
~ -0.8pi + 0.6p4 (52) 

with ~ 0.02. Clearly, due to the numerical nature 
of the Fisher matrix calculation, some sub-dominant con¬ 
tributions in the directions of p2, Ps and persist; we 
ignore these. We see here that o-i, the better constrained 
of the two parameter space directions, is essentially a 
weighted sum of the directions pi and p 4 . di is the vector 
orthogonal to di, and consists of a weighted difference. 

In order to interpret the level of constraint placed on 
the eigenvector directions, we consider equation |49l Re¬ 
call that because we have chosen an ansatz for 5pi(x) 
which is non-negligible at late times only fequation I50L 
E(a;, k) is the dominant non-GR-|-ACDM contribution to 
while only fl{x,k) contributes to fcr^. In exam¬ 
ining equation 1491 it is clear why a sum of directions pi 
and p 4 is forecast to be better constrained than a differ¬ 
ence: both E(x, k) and p,{x, k) are directly sensitive to 
weighted sums of Spi(x) and Sp 4 (x). 

We now consider eigenvectors in the subspace of P 2 , P3 
and P5. These are given by: 

a3 = O.OOOSpi -|- 0.3p2 — 0.8p3 + 0.0009p4 + O.Sps 

~ 0.3p2 — 0.8p3 + O.Bps (53) 

with ~ I0“®, 

d4 = —O.OOdpi — 0.8p2 + 0.1p3 — 0.002p4 -I- O.Gps 

~ —0.8p2 + O.Ip3-I-0.6p5 (54) 

with ~ 4 X 10“"^, and 

ds = —1 X 10“®pi -I- 0.6p2 + 0.6p3 -I- 5 X 10~^p4 -I- O.Gps 
~ 0.6p2 + 0.6p3 -I- O.Gps (55) 


with vxT^ ~ 2. Recall that it is not meaningful to directly 
compare these values with those for eigenvectors di 
and a2. 

In order to understand the relative constraints on the 
eigenvectors in this subspace, we first find the relation¬ 
ships between parameter values along the eigenvector di¬ 
rections. We then use these relationships in equation 0^ 
to find simpler expressions for E(x, fe) and fl(x^k) : 


In the direction d'3, p® — ”^4*3 P 5 — “fPs- 

Therefore, E{x,k) ~ and 

FjIt k) ~ q Ppd 

P{X,K) — O(20ffo)^ n««(a:=0) ' 


• In the direction a4, P 2 — —Sp® and pg ~ Gpg. 

- -i (2ogoVnpffo) ^"’ and p(x,k) ~ 


,pld ,2 


• In the direction cJs, Pq — P3 and pg — pg. Tj{x, k) ~ 
0, and /2(x, k) ~ 0. 


With this manipulation, it can be seen why the direc¬ 
tion d 3 is better constrained than d 4 . The expressions 
for S(a;, k) demonstrate that a small change in param¬ 
eter values along direction as has a numerically larger 
effect on the value of Y,(x,k) than does a small change 
along direction cf4, by a factor of Although p(x, k) is 
affected slightly more by a change along the a4 direction 
than along the ds direction (by a factor of |), this effect 
is clearly subdominant. It is also plain to see why the 
direction dg is by far the worst constrained of this set. 
As long as P2 = P3 = Pg, both I](x, k) and pl ( x , k) are en¬ 
tirely insensitive to the value taken by these parameters. 

It is tempting to attempt a comparison of the above 
results directly with those of [IJ. However, this would 
be misleading, as the Principle Component Analysis 
employed in that work allows the functions Pi(x) to 
take any time-dependence, whereas we restrict the time- 
dependence to that given in equation 1501 Additionally, 
the authors of choose to impose a prior upon the vari¬ 
ance of each function Pi(x), enabling them to meaning¬ 
fully discuss constraints on directions which combine the 
two parameter subspaces which we have considered. Nev¬ 
ertheless, there is some small comparison we can draw: 
in it was found that the functions Pi{x) can not be 
individually constrained, and we similarly find no evi¬ 
dence of any well-constrained direction corresponding to 
a single Pq. 


V. CONCLUSIONS 

The goal of this work has been to understand how weak 
leasing measurements are affected by the individual phys¬ 
ical effects of alternative theories of gravity. In this spirit, 
we have chosen to prioritise clarity. We have therefore 
restricted ourselves to considering gravitational param¬ 
eters, and have not included uncertainty in galaxy bias 
or in the standard cosmological parameter values at this 
time. 

We have constructed an expression for the power spec¬ 
trum of the weak leasing observable convergence, as given 
in eauation l29l By considering only small deviations from 
GR-I-ACDM, we have derived an expression which sep¬ 
arates into an integral over two terms: a general rela¬ 
tivistic kernel, and a source term which encompasses all 
deviations from general relativity. This source term is 
composed of additive terms which are themselves each 
representative of a different effect due to modifying gen¬ 
eral relativity. This neat separation, first of the full ex¬ 
pression into kernel and source and then of the source 
expression into physical effects, allows degeneracies be¬ 
tween gravitational parameters to be physically inter¬ 
preted. We note also that our expression is reliable in 
the ACDM-like case of /3(x) = 0, whereas complimen¬ 
tary works employing a fluid model have found this limit 
difficult to constrain (d^ . 

With equation[23in hand, we first investigated degen¬ 
eracies in the simplified case of a ACDM-like expansion 
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history, showing how the degeneracy direction of weak 
lensing in the /2o — Sq plane relies on the time-dependent 
ansatz for /2(a;). However, we also demonstrated that 
even for an extreme ansatz, weak lensing and redshift- 
space distortion measurements remain non-degenerate in 
this plane, and hence are a viable combination of observ¬ 
ables to offer joint constrains on fiQ and Eq- 

We then moved on to exploit the potential of our ex¬ 
pression as a valuable tool in conducting and interpret¬ 
ing Fisher forecasts of weak lensing observations. We 
found that the linearity of our expression in the gravi¬ 
tational parameters meant that it was technically sim¬ 
ple to compute the required Fisher matrices. Perhaps 
more importantly, the simple form of our expression was 
also found to provide physical interpretation of the de¬ 
generacies which presented themselves in the forecast of 
multi-dimensional constraints. 

We first demonstrated this use of our expression by al¬ 
lowing the effective equation of state of the dark energy 
component to take a CPL form, given in Section llVBI We 
found that in the case where Wa = 0 and a phenomeno¬ 
logical ansatz for fi(x) and E(a;) is assumed, forecast con¬ 
straints on /2o from weak lensing and redshift-space dis¬ 
tortions are highly sensitive to the level of constraint on 
wq from BAO. Forecast constraints on Eq, on the other 
hand, were found to be nearly unaffected. The separation 
of our expression into source and kernel made evident the 
stark difference between the redshift dependences of the 
kernels of and fas- This provided a clear expla¬ 
nation for the higher sensitivity of redshift-space distor¬ 
tions to early time effects, and hence to changes in wo- 
Relaxing the requirement that Wa = 0, we then saw that 
constraints in the p-o — Eg plane were only moderately af¬ 
fected by marginalising over Wa- This offers the exciting 
prospect that the analysis of data from a DETF4-type 
survey may be able to allow for a time-dependent back¬ 
ground expansion without sacrificing much constraining 
power in the po — Eg plane. 

Finally, we used our expression for in combination 
with that for fas , to understand the best-constrained di¬ 
rections in the space of parameters of the scale-dependent 
ansatz given in [l^ . By varying the parameter values 
along the eigenvectors of the Fisher matrix within the 
source terms of and of fas, we could understand an¬ 
alytically why different directions in the parameter space 
were forecast to be well- or poorly-constrained. Using 
again a linear response approach, we were able to explain 
why a weighted sum of the parameters not associated 
with scale-dependence is forecast to be well-constrained, 
and why a direct sum of the parameters associated with 
scale-dependence is totally unconstrained. This work ex¬ 
pands upon the work of [2^ , both by forecasting for a dif¬ 
ferent survey, but also by providing clear reasons for the 
forecast constraints found, something that is non-trivial 
to do using a Principle Component Analysis method. 

The methods we have developed to better understand 
the effect of altering the theory of gravity could in prin¬ 
ciple be extended to understanding degeneracies in other 


scenarios of observational cosmology. We hope that we 
have provided here the groundwork for such develop¬ 
ments. As previously noted, in the interest of clarity 
we have not included uncertainty on the galaxy bias or 
on standard cosmological parameters in our forecasting. 
Ideally, these would be included to provide more accu¬ 
rate forecasting and more general conclusions. However, 
in seeking to include these additional parameters, the 
clarity of equation [29] is compromised, and as such we an¬ 
ticipate that future work in this direction will necessarily 
involve the sacrifice of the descriptive power shown here. 
In such future work, forecasting methods involving sam¬ 
pling the posterior of the joint probability distribution of 
the relevant parameters are likely to be more appropriate. 
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Appendix A: Derivation of 'H^UmImg = 


Here we derive a) the formula for ST-Lix) used in equa- 
tion|23l and b) the relation = 'H"^^m\gr used 

in Section UTOl 

First we write the Friedmann equation as: 


[pM{a) + PD{a)\ 


3 

SttG 


PMoe 


' PDoe 


— J dx' {l+3wo{P)) 


3 


(Al) 


where in the second line we have changed the independent 
variable to a; = In a, and used the energy density evolu¬ 
tion for a fluid with a general equation of state ujd- In the 
third line we have defined R = flAo/flMo, identifying the 
present fractional energy density in the dark fluid (floo) 
with that of the apparent cosmological constant today 
(Had)- We have also made use oi wd{x) = — l + (i{x), as 
introduced in equation 1131 

We define u{x) = J^P^')dx' as in equation [HI and 
assume as in Section III Bl that for a viable cosmology 
|,5(a;)| 1 and |M(a:)| -C 1. Then, expanding the expo- 
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nential and taking a square root leads to: 


n{x) = Ho^/^e-^ [1 + Re^^il - 3u(x))] ’ 

3 f u{x)Re^^\ 


i^o\A^e-t [1 + Re^^y 


1 - 


= R-GRix) 


1 - 


3 f u{x)Re^^'\ 
Re^^ ) 


1 + 


1 + Re^^ ) 
(A2) 


where in the second line we have used a Taylor expansion 
to linear order. Finally, we use the (easily derived) result 
that in ACDM, and assuming negligible radiation, 


(1 - = 


to obtain: 


Re^^ 

1 + Re^^ 


(A3) 


dnix) = nix) - RGRix) = -RHgr{x)u{x) (l - . 

(A4) 


and slip relation, but they may generally still have a dif¬ 
ferent value than 'i!Qfi(x,k) and ^GR{x,k), due to any 
difference in the history of the growth of overdensities. 
Now, from equation [TH] above, we see that 

<f>(a:, fc) = (1 + 5pl{x, k)){l + SAix, k))^GR{x, k). (B2) 

In fact, $Q^(x, fe) = (1 + 5a{x, k))^GR{x, fc), so that we 
have: 


^{x,k) = {1 + 5^{x,k))^Gji(x,k). (B3) 

We also know from equation [3] that in the quasistatic 
regime, 3>(x, fe) = (1 + 5^{x, fe))'h(x, fc). Therefore: 


'I'(a:, fe) = 


1 -I- 5/x(a:, fe) 
1 -I- 5^{x, fc) 
1 -I- 5/x(a:, fc) 


$Q^(a:, fc) 


l + (57(x,fc) 

~ (1 -k 5^l{x, fe) - 5-f[x, k))'l!%ii{x, fe) 


(B4) 


With this result in hand, it is simple to show that the 
combination does not change under the perturba¬ 

tions about the GR+ACDM model. Expanding to first 
order about the fiducial model: 

= [Rgr + <5'H)^ + S^Im) (A5) 

^ n^R^fP (i + 2^ + + om^). 

(A6) 

We call upon a result derived in Appendix A of dl (for 
brevity’s sake we will not repeat the derivation here): 

SUm = Mx)^m^ [1 - ■ (A7) 

Substituting equations (HU and (IA7|) into equation (IA6I) 
one finds 'H'^VLm = RPbPPPi as stated in the text. 

Appendix B: Converting between {5ii{x,k),5^{x,k)} 
and {ji{x,k),Y,{x,k)} 


where the second line comes from the fact that 
$g^(x, fc) = '^Qj^{x,k). Referring to equation IBli we 
see that 


fi{x, fc) = S^j,{x, fc) — <57(0;, fc). 


(B5) 


as in equation 1371 

From here, we can use equations IB3I and IB4I to write 
(suppressing time- and scale-dependence for brevity): 


'k -k $ = 


1 +'^fc‘.T,S 


1 -k <57 


^GR + (1 + ^k‘)^GR 




1 -k <57 


^GR + (1 + ^k-)^GR 


1 + Sfj, 

1 -k i57 


-k (1 -k S^) 


1 


+ ’^gr) 


~ ( 1 -k (5/r - -(57 ) {^gr + '^gr) 


(B6) 


and by comparison with equation IBli we have that 


We derive here the relationship, given in equation 
[371 between two sets of functions which parameterise 
deviations from GR-kACDM in the quasistatic limit: 
{S^(x,k),S^{x,k)} and {p,{x,k),Tj{x,k)}. In doing so, 
we will refer heavily to equations 4 and 5 of For ref¬ 
erence, we reproduce them here, changing to our time 
variable x = ln(a) and making some minor notational 
alterations in order to maintain our conventions: 

^'(a::, fc) = [1 -k p,{x, fc)]4'gfl(x, fc) 

^'(a::, fc) -k $(a:, fc) = [1 -k E(a:, fc)](4'g^(a::, fc) -k fc))- 

(Bl) 

We have labeled the general relativistic potentials with 
an S due to the fact that they are slightly different from 
what we call '^GR{x,k) and 4>G_R(a:, fc). 4'Q^(a:, fc) and 
^QR{x,k) follow a general relativistic Poisson equation 


T,{x,k) = S^{x,k) — -S'-f{x,k) (B7) 

as in equation 1371 

Appendix C: Derivatives of and fas 

We present here the derivatives of PR and fas which 
are required for forecasting in Section llVI These are de¬ 
termined in a straightforward manner by differentiating 
equations [29l and [361 For brevity, we will use IC{x,i) as 
defined in equation |3T] to denote the general relativistic 
kernel in the definition of PR{i), and Gf{x, x) as defined 
in equation 34 of [1^ to denote the general relativistic 
kernel in the case of redshift-space distortions. 

First, the derivatives with respect to fcto Eq are 
given as: 
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5/2o 

d^o 

dfcrs{x) 

dflo 

df(7s{x) 

5So 


dx IC{x, i) 


• —OO 

rO 


= 2 / dxlC{x,i) 


3 / dxVt^{x)I{x,x) 

. ^ — OO 

nf^ix) 


nf^{x = o) 


nf^{x = o) 


Gf{x,x) 


fl^^{x = 0) 


dx 
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(Cl) 


Next, derivatives with respect to wq and Wa- 
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I 3 /■“ 3t 
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Finally, derivatives with respect to {pq}: 
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dfas{x,k) 

dPo 


[ Gf{x,x') 

*/ —OO 


n<l^{x' = 0 ) 


k^dx'. 
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